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Abstract 

Network-assisted device-to-device communication is a promising technology for improving the 
performance of proximity-based services. This paper demonstrates how the integration of device-to- 
device communications and dynamic time-division duplex can improve the energy efficiency of future 
cellular networks, leading to a greener system operation and a prolonged battery lifetime of mobile 
devices. We jointly optimize the mode selection, transmission period and power allocation to minimize 
the energy consumption (from both a system and a device perspective) while satisfying a certain 
rate requirement. The radio resource management problems are formulated as mixed-integer nonlinear 
programming problems. Although they are known to be NP-hard in general, we exploit the problem 
structure to design efficient algorithms that optimally solve several problem cases. For the remaining 
cases, a heuristic algorithm that computes near-optimal solutions while respecting practical constraints on 
execution times and signaling overhead is also proposed. Simulation results confirm that the combination 
of device-to-device and fiexible time-division-duplex technologies can significantly enhance spectrum- 
and energy-efficiency of next generation cellular systems. 

Index Terms 

D2D communication, energy-efficient network, dynamic TDD, mode selection, power control. 

I. Introduction 

During the last decade, wireless communications have experienced an explosive growth, 
both in the number of mobile subscribers and in the data traffic demands. This phenomenon 
is expected to continue in the future Q, mainly driven by an increasingly rich web content, 
video- and audio-streaming and file sharing. Meeting this increasing demand causes the energy 
consumption of wireless systems to escalate. The growing energy bills of operators, limited 
battery lifetime of mobile devices and environmental concerns are progressively steering the 
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research community towards the development of energy-efficient wireless communications Q. 
Efficiency and scalability are therefore becoming the key criteria for the development of the next 
generation (5G) systems, where Device-to-Device communication (D2D) is recognized as one 
of its promising technology components P|4] . 

A natural question in the context of D2D communication is under which condition two users 
should communicate through a direct link rather than via the BS. This problem is known as the 
mode selection problem. The optimal mode selection depends on the performance measure that 
we wish to optimize. For example, the authors in Q select the communication mode to maximize 
user rate while satisfying SINK constraints on active cellular links. In Q, the authors focus on 
maximizing the power-efficiency of the network. Maximizing the quality-of-service (QoS) for 
a given transmit power (and the dual problem of minimizing the power for a given QoS) was 
considered in 0 . To realize the full potential of D2D communications, the mode selection 
should be done jointly with the radio resource allocation. Several works have investigated this 
joint problem [[8]|-[[T^, mainly developing mixed-integer programming models and dealing with 
NP-hard formulations. Hence, heuristics pT||-p^ and game-theoretic approaches plj | have also 
been explored as more practical alternatives. 

Typically, state of the art literature deals with minimizing the transmission power consumption 
or maximizing the overall system throughput. In this paper, we focus on the energy consumption. 
Specifically, we analyze the energy reduction that can be obtained by combining D2D communi¬ 
cation with dynamic Time-Division Duplex (TDD) technology, which is considered an attractive 
duplexing method for 5G networks | |15|16] . In dynamic TDD systems the BS can adjust the time 
allocated to uplink (UL) and downlink (DL) traffic dynamically. Traditionally, this feature has 
been used to compensate for the asymmetry in the UL and DL traffic demands. In contrast, we 
use this degree of freedom to minimize the energy consumption of the transmissions, leveraging 
on the observation that the energy required for sending a fixed amount of data decreases with 
the transmission duration (see [ [T7] | and references therein). 

To the best of our knowledge, the energy efficiency improvements that can be obtained by 
integrating D2D and flexible TDD has not been investigated in literature. Frameworks for D2D 
enhanced TDD networks are proposed in [[T8]|-[[20||. However, they do not account for the mode 
selection and they mainly focus on the adaptive UL/DL slot allocation to D2D pairs, so to 
balance the traffic load, coordinate the interference, improve coverage probability and sum-rate. 
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This paper extends our previous study [ |2T] | to a more realistic multi-link scenario. We consider 
two possible spectrum allocation strategies: one where each D2D pair is assigned its own 
frequency channel, and one where all D2D connections share the same channel. We develop 
algorithms that minimize the energy consumption from the perspective of both the entire system 
and of the mobile devices, since these are the most energy-sensitive part of the network. 

Our main contributions can be summarized as follows: 

• We formulate the joint mode selection and resource allocation problems as a mixed-integer 
nonlinear programs (MINLP), which are NP-hard to solve, in general. 

• For the interference-free case, we propose a low-complexity algorithm that exploits problem 
structure and finds the optimal solution in polynomial (and in some cases linear) time. 

• When considering interference, finding the optimal solution becomes much more challeng¬ 
ing. However, we design an algorithm based on branch and bound (B&B), which computes 
the optimal solution in a much more efficient way than a naive exhaustive search. 

• We propose a heuristic algorithm for computing near-optimal solutions while respecting 
practical constraints in terms of execution times and signalling overhead. 

We structure the paper as follows. § |I^ describes the system model and states some basic 
assumptions. §[115 presents the general problem formulation. §|IV] elaborates the optimal solution 
for the interference-free case, while § |V] presents both the B&B and the heuristic approach for 
minimizing total device energy when all D2D communications share a single frequency channel. 
Numerical results are presented and discussed in § |V^ Finally, § |VII| concludes the paper. 


II. System Model and Assumption 

We consider a single-cell network where infrastructure-assisted device-to-device communica¬ 
tion is enabled. Communication between in-cell users is done in one of two possible modes: 

1) Cellular mode: the transmitter first sends the data to the BS, which then forwards the 


message to the intended receiver (Fig. |l(a)| ); 

2) D2D mode: a dedicated direct link between the transmitter and receiver is set up (Fig. |l(b)T ). 
In the cell, we consider a set C of user pairs that wish to communicate. Each user pair 
constitutes a logical link that we label by an integer 1,2,... L. The BS is denoted by 0 and we 
refer to the users in pair-/ as transmitter-/ (Tx-/) and receiver-/ (Rx-/), respectively. 
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Fig. 1. Communication modes: user pair-/ can communicate either via the BS (a) or through a direct link (b). 


The system bandwidth is divided into orthogonal channels of size W Hz, and time is divided 
into frames of fixed length of T seconds, see Fig. The BS manages the spectrum by assigning 
a time-frequency physical resource block to each logical link. Each resource block consists of 
one frame in the time domain, and one channel in the frequency domain. 
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Fig. 2. Frequency-time resources configuration for D2D communications with dynamic TDD scheme. D2D links can use the 
full frame duration T; UL and DL durations for cellular communications can be reconfigured at each time frame. 


We consider overlay in-band D2D communication, where D2D communications are allowed 
to use cellular resources but may not cause interference to traditional cellular communications. 
To this end, we assume that the channels are divided into two non-overlapping sets: one allocated 
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to the pairs communicating in cellular mode and one to those communicating in D2D mode. 


A. Communication in cellular mode 


To prevent intra-cell interference between concurrent transmissions, the BS follows the channel 
allocation policy of legacy LTE systems and assigns a separate channel to each user pair in 
cellular mode. We assume that the system adopts a dynamic TDD mode, where the UL and 
DL transmissions for a user pair occur on the same frequency channel but alternate in time 
(see Fig. |^. In particular, we consider a dynamic TDD system, where the time allocation 
for UL and DL transmissions can be reconfigured in each time frame, but it is the same for 
all communications. This synchronization of cellular communications is a standard practice in 
current TDD systems to avoid strong inter-cell interference P2|23| |. We denote by fui 
the portion of the time frame allocated to the UL and DL transmissions, respectively. 

Let pio and poi be the transmit power levels used by Tx-l in the uplink to the BS, and used 
by the BS in downlink to Rx-l, which are subject to bounds and respectively. The 
instantaneous rates r/o and tqi achieved in UL and DL, respectively, are assumed to follow 
Shannon’s capacity formula: 


rio = W log + roi = Wlog(^l + ^^y (1) 

Here, Gio is the channel gain between Tx-l and the BS, Gqi is the channel gain between the 
BS and Rx-l, and cr^ is the noise power. The maximum instantaneous rates (corresponding to 
maximum power transmissions) in ([T]) are denoted by and respectively. 


B. Communication in D2D mode 

In D2D mode, each pair can use the full frame duration for its single-hop (p^) transmission, 
as illustrated in Fig. Let U <T denote the active time of pair-/ in D2D mode. As discussed 
above, the D2D communication occurs on different orthogonal frequency channels than those 
used for cellular communications. We consider two channel allocation strategies for D2D pairs: 

• Full Orthogonality (FO): All D2D communications are assigned orthogonal frequency 
channels. Hence, no receiver is interfered by other transmissions within the cell. 

• D2D Resource Sharing (RS): All D2D communications are assigned to the same frequency 
channel, hence interfere with each other. 
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Let Til denote the instantaneous rate between users of pair-/ when transmitting in D2D mode: 

ra = W\og(l + ^y ( 2 ) 

Here, Gu is the direct channel gain between Tx-l and its intended receiver, pu is the transmit 
power level (upper bounded by and Ii = Yhj^iVijGji is the interference power level 

experienced at Rx-l, due to other concurrent D2D transmissions. Note that // = 0 if we consider 

the FO channel allocation scheme. We use to denote the maximum instantaneous achievable 
rate corresponding to 

C. Rate constraint, power feasibility and energy cost 

To guarantee a certain QoS, each pair-/ has a traffic requirement of bi nats during a time 
frame, irrespectively of the communication mode. This QoS requirement can be translated into 
a session rate requirement of hfT nats per second. Specifically, if pair-/ is in cellular mode, 
the transmission times for UL and DL, along with the corresponding instantaneous transmission 
rates, must satisfy 

fiolui^bi and roitdi>bi, (3) 

Similarly, if pair-/ is in D2D mode, f and ru must satisfy 

niU > k. (4) 

The limitation on the transmission power levels, together with the session rate requirements 
above, entail the need to verify under which conditions the communication of a pair can be 
supported by the network. To this end, we introduce the concept of power-feasibility: 

Definition 1 (Power feasibility). We say that user pair-l is power-feasible 

(a) in D2D mode if > bp 

(b) in cellular mode if there exists a time allocation (/uh/di) such that 

l-al + /^dl 

< fui > (5) 

WO 

f > 

^dl _ .^max • 

V ^01 

Assumption 1. Every user pair-l in the cell can be supported at least in cellular mode. 
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Remark 1. The power-feasibility condition Q implies that must satisfy 


^max — “ — ^max ' 

WO ' 0 / 


( 6 ) 


Thus, As sumption^of a feasible UL time allocation for all pairs suggests the following feasibility 
condition for the session rate requirements: maxl^l^r} < min{T — 

By inverting the power-rate relationships ([T]) and we find the energy required to satisfy 
the session rate requirement bi on link-/ for a given time alloeation 

UL, eellular mode; 

DL, cellular mode; (7) 
D2D mode. 


Eioitui) = Piotui = (^exp 

Eoi{tdi) = poddi = ^exp 
=PiiU = fexp 


wu 

bi 

Wtai 

k 


1 ) ^ /^ul) 


a 


- 1 ) y^tdh 
^01 

- 1 1 


WtiJ J Gu 

These functions are convex and monotonically decreasing p7] | in their arguments. This obser¬ 
vation leads to the following result: 

Lemma 1. Any energy-optimal solution must allocate the full frame duration for communication. 
For the D2D mode, this implies that 


min Ell (//,//) = min (T, //), 

H,h H 

while for cellular communication it must hold that 

f ul + tA\ — T. 


The energy cost for communicating in cellular mode includes the energy cost of both the 
transmitting device and that of the BS. However, the BS often has access to cheap and abundant 
energy in comparison with the user equipment, in which case it is relevant to only focus on the 
device energy. To this end, we consider the following two definitions of energy consumption for 
a generic user pair-/ in cellular mode: 

• The System Energy consumption (SE): is the energy consumed by both Tx-/ in UL and 
the BS in DL. By Lemma [TJ the total energy cost is obtained by minimizing 

Eu^^^itui) = Eiofui) + Eoi(T — /ui), (8) 

which is a convex function of t^h 










The User Energy consumption (UE): is the energy consumed by Tx-l in UL transmission, 
disregarding the energy spent by the BS, that is 


which is a convex and monotonically decreasing function in tuh 


(9) 


III. Problem Statement 

We consider the problem of jointly optimizing the communication mode, the UL/DL transmis¬ 
sion times and the powers allocated to all transmitters, to minimize the energy consumption of 
the system. For communications in D2D mode, we consider the two possible spectrum sharing 


strategies introduced in § II-B: FO and RS. The energy cost for cellular communications is either 


the user energy (UE), or the total system energy (SE) presented in § |II-C[ Hence, we consider 
the following four variations of the energy minimization problem: 

• Fully orthogonal channels - system energy minimization (FO-SE), 

• Fully orthogonal channels - user energy minimization (FO-UE), 

• D2D resource sharing - system energy minimization (RS-SE), 

• D2D resource sharing - user energy minimization (RS-UE). 

In each case, we develop algorithms that find the optimal mode selection and resource allocation 
that minimize the chosen energy cost. Specifically, the mode selection policy divides the set of 
user pairs C into two subsets: T>, representing the pairs that should communicate in D2D mode, 
and C = C\V, representing the user pairs that should communicate in cellular mode. The 
mode selection is jointly optimized with the allocation of the optimal transmission power to 
each transmitter, and the optimal UL/DL time allocation for all user pairs in C. 

As will be seen in the subsequent sections, the optimal joint mode selection and resource 
allocation problem for each of these four cases can be formulated as a nonlinear mixed-integer 


mathematical program. In § |IV} we show that optimal solution for FO-SE and EO-UE can 
be found in polynomial time, even though the overall problem is not convex. In § |V| on the 
other hand, we will show that finding the optimal solution in the RS-SE and RS-UE cases is 
challenging, mainly due to the large amount of runtime needed and knowledge of all cross¬ 
link gains required. Thus, we further propose a heuristic algorithm that efficiently computes 
near-optimal solutions and takes into account practical implementation aspects. 





9 


IV. Minimum-Energy Mode Selection with Full Orthogonality 

In this section, we show how the jointly optimal mode selection and resource allocation with 
full orthogonality can be found in polynomial time. For ease of exposition, we first derive the 
optimal solution for a single user pair, and then extend the results to multiple user pairs. 


A. Single user pair 

We first characterize the minimum energy cost for cellular and D2D communication. We show 
that the optimal system energy (SE) cost can be found by solving a simple convex optimization 
problem, while the minimal user energy (UE) cost admits an explicit expression. 

1) Minimum energy cost for communication in cellular mode: In cellular mode, the UL/DL 
time allocation is chosen to minimize one of the following two objectives: 

• Minimizing SE: The minimum amount of energy of pair-/ in cellular mode can be deter¬ 
mined by solving the following single-variable convex optimization problem: 


minimize Eiq if^x) + Eqi {T - Ex) 
bi bi 


tul 

subject to 


< fui < T - 


^max — — ^max ^ 

WO '0/ 


(lOa) 

(10b) 


where constraint ( |10b| ) ensures power feasibility in the sense of Definition Problem ( fT0| ) 
can be solved efficiently using a wide variety of methods, including simple bisection 


search p5| . Let denote the optimal solution to ( fT0| ). The energy cost of pair-/ in cellular 
mode is then 


= £,„(<;,) + E„(T - 1;,). d i) 

• Minimizing UE: Here, the only difference from the problem formulation in ( [T^ is that 
the objective function reduces to Eio{Ex). By monotonicity of the objective function, the 
optimal solution is attained by /*) = T — with the corresponding optimal energy cost 

= E,„{T - (12) 

'01 

2) Minimum energy for communication in D2D mode: In the D2D mode, no traffic is for¬ 
warded through the BS and only the user equipment consumes energy for the connection. The 
minimum energy cost follows from Lemma [TJ with f = 0: 

- (-» {m) - 0 
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Equation ( fT^ is only valid when the D2D mode is power feasible. Since we have ensured power 
feasibility only for communications in cellular mode, we need to verify that r™^T < bi before 
applying ( fT^ . It turns out to be more convenient to work with the extended value function 


pD2D 
^ll = 


(wt) 0 


2 

rp 


Gu- 


(14) 


if < bi 

+00 otherwise. 

3) Optimal mode selection policy and resource allocation: The optimal mode selection policy 
consists in first solving the convex optimization problem to estimate the energy cost for cellular 
mode, and then comparing it with the energy cost for D2D mode. In the FO-SE scenario, we 
compare ( fT4| ) with ( [TT] ), and in the FO-UE scenario we compare ( fT4| ) with ( fT2| ). The optimal 
communication mode is simply the one that requires the least amount of energy. 

Once the optimal communication mode and transmission time has been selected, the corre¬ 
sponding optimal powers are easily derived as 

hi 


pio = exp 


wt:, 


- 1 


a 

Gio 


Poi = exp 


WiT - 


a 

Gqi 




Pii = exp 


hi 

WT 


-11^, 

Grii 


UL, cellular mode; 
DL, cellular mode; 
D2D mode. 


(15) 


It is possible to interpret the mode selection policy in terms of the channel gain (and, so, 
the physical distance) between the two communicating devices. We will explore this geometrical 


interpretation in § VI to characterize regions in the cell where D2D communication is preferable. 


B. Multiple user pairs 

Under full orthogonality, the main challenge with multiple user pairs is that all communi¬ 
cations in cellular mode must use one common UL/DL time allocation. In other words, the 
optimal UL/DL time allocation must account for the energy consumption of all users. Under 
full orthogonality, D2D links do not interfere with each other. Thus, the minimum energy cost 
for each pair in D2D mode is the same constant value ( [T?] ) as in the single user pair case. The 
energy consumption in cellular mode, on the other hand, depends on the UL time allocation. 
Let m G (0,1}^ denote the mode selection vector whose entries satisfy 


mi = 


0 if pair-/ is in cellular mode, 
1 if pair-/ is in D2D mode. 


(16) 
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The general formulation of the joint mode selection and time allocation problem for multiple 
pairs under both FO-SE and FO-UE cases, is the following MINLP: 


L 


minimize 

m,tui 

- mi) 

1 = 1 


(17a) 

subject to 

—^- Tnii < < T -^—h Tmi, 

^max '' — — ^max '' ’ 

WO '0/ 

V/, 

(17b) 


/ui £ [0)^]) mi G {0,1}, 

V/. 

(17c) 


The objective function is the total energy consumption of all the L user pairs, with 
given by @ or @ under FO-SE or FO-UE, respectively. Constraints ( |17b| ) ensure that pair-/ is 
assigned to cellular mode, only if it is power feasible in the sense of Definition [TJ Since the user 
pairs in D2D mode do not interfere with each other, it turns out to be convenient to eliminate 
the integer variables by incorporating the optimal and feasible mode selection solution in the 
single pair case. Thus, problem ( fTT) ) can be recast as 


minimize F{Ei). 


(18) 


Ei{t 


ulj 


Here E{tu\) — Ei{Ei), and E'/(tui) denotes the minimum energy-cost for the single pair-/ 
when the UL time is fixed to that is 

min{£“'>, ‘■‘■(tul)} if <„l £ [=fe, T - 

otherwise. 

Equation ( fT9l ) reveals the piecewise nature of Ei{tu\). For < bi/r^^^ and E\ > T — 
we have Ei{E\) — Eji"^^, which is a finite constant if user pair-/ is power feasible in D2D mode 
and +(X) otherwise. In the interval T — Ei{t-a\) is equal to the constant 

or given by the function depending on whether and at which points the graphs of 

and intersect. Note that the two graphs can intersect only once if 

is monotonically decreasing (UE) or twice if it is convex (SE). See Fig. for an illustration. 

To better describe the piecewise nature of £'/(/ui), we introduce A/ = as the 

interval of Ei during which Ei{tu\) = If, for a pair-/, such an interval does not exist 

(Le., A/ = 0), D2D mode is always more energy efficient than cellular mode for that user pair. 

Figures and show the minimum energy-cost of each pair, obtained as illustrated in Fig. 
and the corresponding E{tu\) of a simple network with three user pairs, under FO-UE or FO¬ 
SE, respectively. Note that function E{tu\) is non-convex on the interval [0,T]. The following 
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_ ^CELL 


_ ^D2D 



L . 

. i 

min 

^max 



CD 


0.2 0.4 0.6 0.8 , 

^ui T - i 

'’’ai 


(a) FO-UE single link. 


(b) FO-SE single link. 


Fig. 3. Deriving under both the FO-UE case (a) and the FO-SE case (b). The interval depends 

on whether and at which points the graphs of and E^‘^^ intersect, but also on the interval [-flW,?"-mW] for 

power feasibility in cellular mode. The two intervals coincide when the graphs intersect outside the power feasible interval (a). 

A^ C -if at least one of the two possible intersection points is within the power feasible interval (b). A^ = 0 

'^10 ^01 

if the two graphs never intersect. 


lemma establishes the key property of which will be useful later when solving the problem 

formulated in ( fT^ . 

Lemma 2.(ci) In the FO-SE scenario, F(F\) is a piecewise convex function. 

(b) In the FO-UE scenario, F{F\) is a piecewise decreasing function. 

Proof: For each pair-/, by its definition we have that if A; = 0, is a constant 

value on [0,T], otherwise, is a constant value or -foo in the two intervals [0,t™™) 

and (r™'^^,T]. During the interval in the FO-SE case is given by @, 

which a convex function of F\, while in the FO-UE case, is given by 0. which 

is a monotonicahy decreasing function in Fh The function F{F\) is obtained as the sum of 
of all L pairs. Hence, the whole interval [0,T] is divided into J <2L -\-l adjacent 
intervals. In the FO-UE case, E{tu\) is the sum of constants and convex functions in each 
interval, which makes F(fui) piecewise convex. In the FO-UE case, on the other hand, F(fui) 
is the sum of constants and monotonicahy decreasing functions in each interval, which makes 
F(tui) piecewise decreasing. ■ 
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(a) FO-UE three links: energy-cost functions. (b) FO-UE three links: F{tu\). 


Fig. 4. FO-UE problem: minimum energy-cost functions of three transmitter-receiver pairs (a), and their corresponding sum 

F{Ui) (b). 




(a) FO-SE three links: energy-cost functions. 


(b) FO-SE three links: F{tui). 


Fig. 5. FO-SE problem: minimum energy-cost functions of three transmitter-receiver pairs (a), and their corresponding sum 
F(iui) (b). 


Based on Lemma the optimal solution to ( fTSj ) can be computed efficiently. 


Proposition 1. a) In the FO-SE case, let be a partition of [0,T] induced by the points 

{r™", , r™™, Then, the optimal UL time allocation f*) can be found by solving 
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at most 2L — 1 single-variable convex optimization problems of the form 



b) In the FO-UE case, G Moreover, i/max{T;™} < 


then = min{r - 

l ^ ^ l ^ ^ 0 ! 


.max 


Proof: Since for each user pair-/, reaches its maximum value during the two 

intervals [0,rp“) and (T™'^^,r], thus F{tu\) achieves its maximum value in the two intervals 
Fi = [O, mini{T™™}] and Fj = [maxi{r™'^^}, T]. Hence, is not in Fi and Fj, but must be 
found in one of the remaining (at most) 2L — 1 intervals. By Lemma in the FO-SE case, 
F{tu\) is piecewise convex. Hence, its global minimum can be found among its 2L — 1 local 
minima in each interval. In the FO-UE case, from Lemma we know that F(/ui) is piecewise 
decreasing. Thus, its global minimum can be found in the set However, for 

each there is at least one component in the sum defining that decreases for 

tu\ > Therefore, the global minimum can only be found in the set Furthermore, 

if max{r™'"} < min{r™'^^}, then = min{T™®^}. ■ 

Given the optimal solution t*j to problem ( fTS] ), the optimal mode selection vector m* of 
Problem ( fTT) ) is then given by setting, V/ G £: 




1 otherwise 


and the corresponding optimal transmission powers are derived as in ( [TS] ). 

V. Minimum-Energy Mode Selection with D2D Resource Sharing 

To increase the cell capacity and the spectral efficiency of the system, we now consider the 
D2D resource sharing (RS) strategy, where all communications in D2D mode are assigned the 
same channel resource. We show that the optimal mode selection and the optimal power/time 
allocation are harder to compute when interference is considered among multiple communica¬ 
tions. Nevertheless, we develop a combinatorial optimization algorithm that is guaranteed to find 
the optimal solution, and often does so very quickly. This optimal method is complemented by 
a fast heuristic, suitable for real-time implementation under practical signalling constraints. 
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A. Optimal resource allocation via mixed-integer nonlinear programming 

The effect of interference on the energy consumption in D2D mode appears explicitly in Q. 
Due to the fixed time allocation T (see Lemma [T]), minimizing the energy consumption of D 2 D 
communications is equivalent to minimizing the transmission powers. To meet the minimum rate 
requirement in the transmission power of any pair-/ in D2D mode must be such that 


Pii > 


exp 


WT 


- 1 


Gn 


Introducing 7 *®* = [exp (^) — l] as the target signal-to-interference-plus-noise ratio (SINR) 

tgt^2 t t G • 

required to satisfy the session rate requirement of pair-/, Pi = and hij = 7 ;® we can 

re-write this inequality as 

Pii >r)i + '^Pjjhij- (20) 

m 

The joint mode selection and power/time allocation problem can now be formulated as the 
following MINLP problem: 


nCELL 

ni 


minimize ^ {Tpii)mi E\ 

—^- Tnii < < T — 


subject to 


(tui)(l - mi) 
hi 


Tmi 


' 10 


01 


di + 


- C{1 - mi) < Pii, 


0<pii< pr^'m, 
m G {0,1}'^, /uiGIOiT], 



( 21 a) 

V/, 

( 21 b) 

V/, 

( 21 c) 

V/ 

( 21 d) 


( 21 e) 


Here, m is the mode selection vector defined in (16). Constraint (21b) ensures power feasibility 


for pairs in cellular mode, while ( | 21 c[ ) guarantees that the rate requirement is satisfied for each 
pair in D2D mode. The constant C in pic| ) is a large number (C = max/{r^/ -|- 
example) ensuring that the constraint is only enforced for users in D2D mode. The expression for 
in the objective function is either ([ 8 ]) or (|^, depending on whether we are interested 
in the RS-SE or RS-UE problem, respectively. 

Problem ( |2T] ) belongs to the class of mixed boolean-convex problem, where for each fixed m G 
{0,1}^ the objective function is convex in the continuous variables. In general, MINLPs are NP- 


hard problems |26|, combining the combinatorial difficulty of optimizing over discrete variable 
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sets, with the challenges of handling nonlinear functions. Its solution time grows exponentially 
with the problem dimension. If L is very small {i.e., smaller than 15), it can be solved exactly 
by exhaustive enumeration of the 2^ possible mode selection vectors. However, realistic cellular 
networks might consist of a large number of user pairs. For this reason, we propose an algorithm, 
based on a B&B strategy (e.g., ]27|28[ ), that often finds the optimal solution to ( |2T] ) in a much 
more efficient way than the exhaustive search. 

Due to space restrictions, in the sequel we will focus on algorithms that solve the RS-UE 
problem, considering that the mobile devices are the most energy-sensitive component of the 
network. However, the same approaches can be applied to the RS-SE case. 


B. A branch-and-bound approach for finding the optimal solution 

We propose to solve the MINLP problem using a B&B approach, where possible mode 
selection vectors m are explored through a binary tree. Each node of the tree (except the root) 
represents a subproblem where one of the mode selection variables mu with / e £, is set to 
either 0 or 1. Each branch corresponds to a subset of the possible mode selection vectors. 

Before proceeding, it is convenient to introduce a definition and two useful propositions on 
the feasibility of the mode selection vector m. Let H be the non-negative matrix with entries 
Hij = hij if I j, and zero otherwise, and the vectors p = {piu'il G G 

£)^ and r) — {rju'il G £)^. For each mode selection vector m, we can define the corresponding 
set of pairs assigned to D2D mode and to cellular mode as and Cm, respectively. Let Am 
denote the L x |I?m| incidence matrix, which is formed by removing the /-th column from the 
L X L identity matrix if mi = 0. We define Hm = A^HAm, Pm = A^p, and 77 m = 

The constraint pib| ) can be written in matrix form as 

(Im - Hm)Pm > and Pm < (22) 


where the inequalities are component-wise, Im = A^Am, and Pm'^^ = A^p™'^^. 

Matrix Hm has strictly positive off-diagonal elements, and we can assume that it is irreducible 
because we do not consider totally isolated groups of pairs that do not interact with each other. Let 


p(Hm) denote the largest real eigenvalue of matrix Hm- From the Perron-Frobenius theorem |29|, 
we have the following well-known proposition: 
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Proposition 2 ( p0| | Chapter 2). For a given mode selection vector m, the necessary and sufficient 
condition for the existence of a positive Pm to solve inequality (Im — H„i)Pm > Vm is that 

p(H^) < 1. (23) 

Moreover, p^ = (1^ — is its component-wise minimum solution. 

Proposition provides an easy condition to verify if a mode selection vector m is feasible. 

Definition 2 (Feasible mode selection vector). A mode selection vector m is feasible if both 
condition ( |2^ and (Im — Hin)“^r7m < Pm^ verified. 

Proposition 3 (pT||). If m. is not feasible, every other mode selection vector m, such that the 
set {I E C : mi = 1} {I E C : fill = 1}> is not feasible. 


The main idea of B&B is to only explore branches of the binary tree that have the potential 
to produce better solutions than the best solution found so far. This is done by computing upper 
and lower bounds on the optimal value at each node. If the lower bound of a node is larger than 
the current upper bound, then there is no need to explore its branches. 

To achieve a good performance of B&B, it is essential to select the branching rule and tree 
exploration strategies carefully, and to have efficient methods for computing good upper and 


lower bounds p2| |. In our implementation, we have made the following choices: 

1) Initial upper bound. By assumption, letting all pairs communicate in cellular mode, i.e. 
setting mi = 0 for all I, is always feasible. We therefore take the corresponding energy cost 
as an initial upper bound on the optimal cost, computed by solving 

(24) 


minimize 


lec 


subject to max {——} < < min{T 


leC 


10 


leC 




0 / 


where is given by (|^ in the RS-UE case. 

2) Branching rule. The branching rule selects the next variable to fix at each node of the tree. 
The goal is to identify the branching variable that changes the problem the most, either to 
quickly detect branches that can be cut, or to significantly improve the current solution. 
Our branching strategy for the RS-UE problem is based on first solving the FO-UE problem 


described in § 1^ The corresponding optimal mode selection vector reveals the 
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set = {/ 


m 


FO 


= 1} of pairs that prefer to communicate in D2D mode in an 


interference-free environment. For each link-/ G we define a measure of its strength 
= Sier>po i^iGii/Gii. This measure attempts to account for both the interference that 
pair-/ produces on the shared frequency resource and its own direct gain. The branching rule 
first selects variables in in order of decreasing si, and then considers the remaining 
variables in an arbitrary order. This rule exploits the fact that the optimal solution of the 
RS-UE problem is often close to the FO-UE optimal, with a few differences in the pairs 
in D2D mode, due to interference; pairs that prefer to be in cellular mode in FO will also 
tend to prefer cellular mode in RS. By fixing pairs with high interference strength first, 
we increase the likelihood of finding infeasible solutions quickly. Once an infeasible mode 
selection is found, we can make use of Proposition and discard all branches below the 
current node in the search tree. 

3) Tree exploration strategy. Once a branching variable nii has been selected by the branching 
rule, the tree exploration strategy determines if the child node to investigate next should be 
in cellular or D2D mode (that is, have mi — 0 or 1). We always prefer to try to allocate 
users to D2D mode first (that is, set the branching variable = 1 first). 

4) Upper and lower bounds. Each node of the search tree corresponds to a partial mode 
selection vector with some components set to 1 or 0, while others are still undetermined. 
We define the two sets of pairs corresponding to the determined (fixed) and undetermined 
variables, as and U, respectively. 

When we consider a node, we first verify if the fixed variables form a feasible mode 
selection vector in the sense of Proposition If they do not, then no bounds are computed 
and the node and all branches below it are disregarded. Otherwise, upper and lower bounds 
are obtained as the sum of the minimum energy cost of the pairs in T and, respectively, 
an upper and a lower bound of the energy cost of the pairs in lA. 

To compute an upper bound, we assign all pairs in lA to cellular mode and we solve 


the problem formulation in ( |24l ) where £ is replaced by /Y. To determine a lower bound, 
we compute the minimal energy cost of the unassigned pairs when they operate in full 
orthogonality ii.e., we solve the FO-UE problem over only the unassigned pairs). To 
understand why this is a lower bound, note that the computation is a proper relaxation of the 
MINLP where the UL times of fixed and unassigned pairs are allowed to differ, and when 
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the unassigned links that end up in D2D mode do not suffer interference. Furthermore, we 
strengthen the lower bound by increasing the noise power of pairs in U by the interference 
that the transmitters fixed to D2D communication in T incur on them. 

Once achieved the optimal feasible mode selection vector and transmission time, the corre¬ 
sponding optimal powers are obtained as in ( [T5] ) for the cellular users, and as in Proposition 
for transmitters in D2D mode. 

The B&B algorithm is guaranteed to find the optimal solution, and does so much faster 
than the exhaustive search, as shown in § |VI-B2[ However, for large networks it can still have 
impractical running times. In addition, the optimization formulation assumes that all cross-gains 
between users are known, something that would require significant communication overhead. 
We therefore turn our attention to heuristics that can be run in real-time and do not assume 
centralized knowledge of all the channel gains. 


C. Heuristic approach to achieve a practical sub-optimal solution 

In this section, we present a heuristic algorithm that achieves a near-optimal solution to ( |^ 
in a more practical and scalable way than the B&B approach. Again, we focus on the UE case. 

The key idea of this algorithm is to first determine an initial mode selection vector, together 
with the corresponding power/time allocation, and then improve this solution by means of a 
distributed power control algorithm based only on local measurements. The heuristic algorithm 
is described in Algorithm [TJ with the two main steps in the following: 


1) Initial phase: We adopt the optimal solution to the FO-UE problem in § |IV] as the initial 
solution, denoted by The FO-UE problem is solved by the BS. 

For each pair-/ G (i.e., assigned to D2D mode), the BS also computes the energy 

it would consume if in cellular mode, that is from and broadcasts 

and E^^^^ to each Tx-/ e T>^fo. 

The initial mode selection vector is obtained under the assumption of no interference 
among the D2D pairs. However, under the RS scenario, all the D2D pairs share the same 
channel, thus can be energy inefficient, or even infeasible, due to the interference. 
Therefore, a distributed power control algorithm is then executed by the D2D pairs to find 
a feasible and more energy-efficient solution: 
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Algorithm 1: Heuristic approach for RS-UE minimization 
Input: {i^,GiuGi^,Goue)\/leC 

Output: m*,p* 

1 (m^®, p(m^®)) ^ solution to FO-UE problem; 

2 each I e acquires (tui(m^®)) from the BS; 

3 p*^°^ p(m^‘^), ^ A: = 0; 

4 each I E E^fo computes 7 ^°^; 

5 convergence True; 


6 while convergence do 

7 E- 


8 

9 

10 

11 


for each I E T> (k) do 


'm^ 

^ tgt 


(fc+i) 7,® (fc) 

Pli =^Pii, 


if p\i~^^^ > min then 


m 


(fc+i) 


^ 1) 


E^(k+i) E- E^{k) \ {/}; 


12 

13 

14 

15 


each I E V^(k+i) computes 7/*'^"''^^ 
if 7f^'^ > 7 ?^*, V/ e P^(fc+i) then 

convergence False; 

p*^p(fc+i)^ 


2) Iterative distributed power control for D2D pairs: Using the iterative power control method 
originally proposed by Foschini and Miljanic in p^ , each Tx-l in D2D mode can achieve 
its target SINK 7 *®* by updating its transmit power as follows 


tgt 

(fc+l) 7/ (k) 

Pli =—)Pu\ 
Ti 


(25) 


where is the perceived SINK for pair-/ G V-oiFo in iteration-/^, defined as 7 ^^^ — 


p\'^Gu 


and p*^°) = p(m^®). 


To achieve a feasible mode selection vector and to further reduce the energy cost, some 
links in D2D mode need to switch to cellular mode. Specifically, pair-/ in D2D mode will 
switch to cellular mode if its transmit power level exceeds its maximum limit or if it is 
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more energy efficient for it to communicate in cellular mode, that is. 


^ min J ^ pCELL f'rv.FO'li r.™ 
Pii > i [tui[ra )) , Pi 


(26) 


where we introduce the design parameter 0 > 1. 

Remark 2. The selection of the parameter 9 accounts for the following key aspects: 

a) Mode switches incur additional signaling overhead between mobile devices and the BS 
to coordinate the re-allocation of radio resources. Hence, if 9 > 1 mode switches will 
occur only if they result in a significant energy gain. 

b) When a pair in D2D mode switches to cellular mode, it requires another orthogonal 
frequency channel for its transmission. The parameter 9 can control this trade-off between 
channel reuse and energy consumption. Specifically, a large value of 9 enforces more pairs 
to communicate in D2D mode and reuse the same channel, even if this comes at the cost 
of a higher energy consumption. 

c) Once the final set of pairs assigned to cellular mode has been obtained, the new optimal 

UL transmission time cannot be larger than the initial which means that the 

energy consumption of each pair in cellular mode could increase. Hence, by setting 9 > 1, 
we reserve a margin for the energy increase in cellular mode due to mode switches during 
the distributed power control process. 

During the power update ( |^ , if Tx-/ finds that condition ( |^ is fulfilled, it asks the the BS 
to switch it to cellular mode and to assign it an orthogonal frequency channel. Otherwise, 


it keeps updating its power according to ( |25] ). The BS keeps track of the pairs changing 
communication mode, and updates the mode selection vector. Specifically, if denotes 
the mode selection vector at the k-th iteration of the algorithm, then 


m 


(fc+i) 


0 , 


if = l,and condition ( [26| ) is verified, 
ml', otherwise, 
where the initial mode selection vector 

This power control algorithm converges to the minimum power levels that the user pairs 
remaining in D2D mode need, to fulfill the rate requirement. Once the algorithm converges, 
the BS recomputes the optimal power/time allocation for the user pairs in cellular mode, 
broadcasting this information before the data transmissions take place. 


(k) 




22 


VI. Simulations and discussion 

This section presents simulation results that validate our theoretical findings and evaluate 
our proposed algorithms. All simulations consider a single cell with a BS, equipped with an 
omnidirectional antenna, positioned in the center of the cell. The simulation parameters listed 
in Table |I| are chosen to represent an urban LTE deployment and are used throughout. 

TABLE I 

Simulation parameters oe the system under study 


Parameter 

Value 

Carrier Frequency 

1 GHz 

Cell Radius 

500 m 

Frequency channel bandwidth (W) 

5 MHz 

Noise Power 

-174 dBm/Hz 

Path-loss coefficient (a) 

4 

Path gain at reference distance of Im (Go) 

5.7- 

Max Tx Power for the BS 

40 W 

Max Tx Power for transmitter-/ (p™'"'') 

0.25 W 

Time Frame duration (T) 

1 time unit 


A. Single link analysis: geometrical interpretation of the optimal mode selection policy 


We begin by developing a geometrical interpretation of the optimal mode selection policy for 
the single link case, under the assumption that the channel gains follow a conventional path 
loss model Gij = GoDf^, where Dij is the physical distance between Tx-i and Rx-j, Go is 
the path gain at a reference distance of Im, and a is the path-loss exponent. Recall that D2D 
communication is preferable when 

Let us first study the user-energy objective, in which case and are given 

by ( [T2| ) and ( fTd] ), respectively. Using the path-loss model, we can transform the mode selection 
policy into the following equivalent condition in terms of the distances between the transmitter, 
receiver and BS: 


/ (pbi/WT _ \ 

(iWrk) 

'-V-' 

«(Do() 


Du < 


Dio — b^{Doi)Dio. 


(27) 
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Note that k depends on Dqi, since Dqi affects and thereby = T — bi/r™i^^. Thus, even 
though we are neglecting the energy cost for the DL transmission, Dqi still influences the optimal 
mode selection. To characterize the region where D2D mode is preferable, we fix the position of 
Tx-/ (and therefore Dio). We then vary the position of Rx-l along a circle centred at the BS, thus 
keeping Dqi and k{Dqi) constant. Inequality ( |T7l ) now states that D2D mode is preferable when 
the distance between the transmitter and receiver of pair-( is less than k{Doi)Dio. In other words, 
D2D mode is more energy efficient when Tx-/ is located in the arc defined by the intersection 
of the circle of radius Dqi centered at the BS, and the disc of radius k{Dqi)Dio centered at 
Tx-L The D2D optimal area can be constructed by tracing out these arcs for various distances 
between Rx-l and the BS, see Fig.|^ One can show that k{Doi), and thus the D2D optimal area, 
decreases as Rx-l gets closer to the BS. 

Fig. 1^ illustrates the D2D-optimal area in red and the D2D power-feasible are in light blue. 
To ensure that Assumption [T] is satisfied, we set bi equal to the maximum traffic rate that can 
be supported when Rx-l is on the cell boundary. 

Although ( |27| ) does not formally describe a disc around Tx-l, the D2D-optimal area is close 
to circular. The reason for this is the power imbalance between the user equipment and the BS, 
which makes very small, T and k 1 practically independently of Dqi. 



Fig. 6. Dashed circles centred at the position of Tx-Z have radius k{Doi)Dio from j27l , and represent, for each of the two 
positions of Rx-l, the area within which D2D mode is more energy efficient than cellular mode for the user energy minimization. 


Similar calculations and arguments can be made for the system-energy (SE) objective. In this 
case, is given by 0 and the expression for k(Dqi) gets a bit more involved 


Fig. shows representative results for our simulation scenario. We observe that the D2D-optimal 
area is no longer circular and that the D2D mode is preferable in a large portion of the cell. 
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0 


500 


0 


- 500 ^ 


- 500 ^ 


(a) Tx-l is 250m away from the BS. (b) Tx-l is 450m away from the BS. 


Fig. 7. D2D optimality area when minimizing the mobile user energy consumption. Red area represents the positions of Rx-l 
for which D2D mode is more energy efficient than cellular mode, while light blue disk represents the area within which Tx-l 
can fulfil the rate requirement transmitting in D2D mode with a feasible power level. 




0 


500 


0 


- 500 ^ 


- 500 ^ 


(a) Tx-l is 250m away from the BS. (b) Tx-l is 450m aeay from the BS. 


Fig. 8. D2D optimality area when minimizing the system energy consumption. Red area represents the positions of Rx-l for 
which D2D mode is more energy efficient than cellular mode, while light blue disk represents the area within which Tx-l can 
fulfil the rate requirement transmitting in D2D mode with a feasible power level. 




B. Multiple link analysis 

The simulation experiments for the multi-link case are set up as follows. We generate random 
network topologies with given number of user pairs. Transmitters and receivers are randomly 
placed within the cell area. An example of network with 10 user pairs is given in Fig. where 
red squares represent the transmitters and green circles represent the receivers, transmitter and 
receiver forming a pair are labelled with the same number. We assume large-scale path loss model, 


where gains are computed as described in Section VI-A Without loss of generality, we further 


assume the same max transmission power level for all mobile transmitters, and the same traffic 
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requirement for all pairs, indicated with b. To ensure Assumption 1 1 we set b = ^mix_)_^max T, where 
and are the maximum achievable rate in UL and DL, respectively, when transmitter 
and receiver are both at the cell edge. For a given number of user pairs, we investigate 1000 
random networks and present the averaged results. 



Fig. 9. Network with 10 user pairs randomly placed in a cell of radius 500m. The red squares represent the transmitters and 
the green circles represent the receivers. Transmitter and receiver forming a pair are labelled with the same number. 


1) Energy gain by enabling D2D communications in a fully orthogonal system: To quantify 
the energy savings that can be obtained by exploiting direct communications, we compare the 
energy cost of the optimal FO-UE solution with the total user energy when all pairs are forced to 
communicate in cellular mode. For each of the 1000 random configurations used in our Monte 
Carlo study, we sort the links in order of increasing energy gain. Fig. shows the averaged 
results for networks with 10 and 30 user pairs, respectively. 

We observe that when D2D communication is enabled, each user either consumes the same 
amount of energy as in the traditional cellular mode or reduces its consumption (and sometimes 
significantly so). Indeed, Fig. [T^ demonstrates that, on average, half of the transmitters in the 
cell achieve an energy gain larger than 20%, and one third of the pairs have an energy gain 
larger than 60%. The average energy saving per user is around 40%. 

It is important to notice that the energy gain is not only a consequence of the proximity of 
users, but also stems from the fact that D2D connections can use the full frame. Thus, it is the 
combination of D2D and flexible TDD technologies that contribute to the total energy savings. 
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(a) Networks with 10 user pairs (b) Networks with 30 user pairs 

Fig. 10. Energy gain by enabling D2D communications in a fully orthogonal system. User pairs are sorted in increasing order 
of the energy gain they achieve by perfoeming mode selection, compared with traditional communication via the BS. 


2) Performance evaluation of the B&B algorithm for RS-UE: The difficulty in solving Prob¬ 
lem ( [2T] ) lies mainly in the possible large search space of integer solutions. Table |I^ shows the 
average number of mode selection vectors explored by different strategies before achieving the 
optimal solution. We consider the naive exhaustive enumeration algorithm, where we only elim¬ 
inate infeasible solutions using Proposition]^ and the B&B algorithm described in Section |WB, 
both with the proposed branching rule and with a random selection of the branching variable. 
Averaged results show that the braching rules has a strong effect on the run-time. 

TABLE II 

Avg. number of explored integer solutions 


Algorithm 

10 user pairs 

15 user pairs 

Exhaustive enumeration 

472.975 

7.46 • 10® 

B&B - Random branching rule 

69.472 

251.21 

B&B - Proposed branching rule 

25.57 

54.72 


3) Performance evaluation of the heuristic mode selection algorithm for RS-UE: In this 
section, we evaluate the performance of the heuristic mode selection policy for the RS-UE 
problem. Fig. [U] shows the additional energy cost of the heuristic relative to the optimal solution 
computed using the B&B solver. For networks with 10 user pairs, the heuristic is within 10% 
of the optimal solution for almost all network configurations (Fig. |ll(a)| ). For networks with 
30 user pairs, the heuristic performs slightly worse. This performance degradation is due to the 
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larger degree of freedom in packing D2D links on the same frequency channel. However, it still 
remains smaller than 10% for most of the configurations. 

In both plots in Fig. there are a few rare network configurations where the heuristic 
performs much worse than optimal algorithm. Fig. is an example of such a scenario. Under 
full orthogonality, both pair 9 and pair 4 in Fig. would be assigned to D2D mode. When the 
heuristic initially attempts to assign these links to the same channel, it encounters an infeasible 
configuration due to the high interference that Rx-4 perceives from Tx-9. Therefore, only one 
of the two pairs can be assigned to D2D mode. The optimal decision is to let pair 4 in D2D 
mode, but during the first iterations of the power control in the heuristic, the interference from 
Rx-9 on Tx-4 leads pair-4 to leave the shared channel and switch to cellular mode. 




Additional energy cost of the suboptimal solution (%) Additional energy cost of the suboptimal solution (%) 


(a) Networks with 10 user pairs. 


(b) Networks with 30 user pairs. 


Fig. 11. Performance evaluation of the heuristic mode selection algorithm for RS-UE. Additional energy cost of the heuristic 
relative to the optimal solution achieved with B&B . 


4) Performance comparison of different mode selection policies: We have developed a number 
of optimal and heuristic mode selection policies. The aim of this section is to compare these 
different approaches and to bring additional insight in how they perform. 

Clearly, allowing for resource sharing reduces the number of orthogonal channels required, 
but increases the energy cost due to the interference among D2D links. Moreover, the heuristic 
allows to penalize the use of channel resources by tuning the threshold value 0. Fig. [T^ shows the 
energy-channel performance for the different proposed algorithms in simulations of networks with 
10 and 30 user pairs, respectively. We note that FO-UE represents the energy-optimal solution, 
using the same amount of frequency channels as the cellular mode. The RS scenario, in which all 
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D2D pairs share the same frequency resource, uses fewer channels but a slightly higher energy 
cost due to the interference between D2D users. The small increase in energy consumption 
can be understood by noting that the transmission powers assigned to D2D pairs are generally 
very small, partly because D2D pairs typically have high direct gains (since the transmitter and 
receiver often are in close proximity of each other), and partly since D2D connections can use 
the full frame. These results demonstrate that D2D communications in cellular networks have 
the potential to improve both spectrum and energy efficiency over a traditional cellular solution. 

Finally, we evaluate the performance of the heuristic method for different values of the 
threshold 6. As expected, large values of 6 decrease the number of channels used at the expense 
of a slightly increased energy consumption. For this reason, 9 is an important design parameter 
for finding a suitable trade-off between energy consumption and channel use. 
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(a) Networks with 10 user pairs. 


(b) Networks with 30 user pairs. 


Fig. 12. Performance evaluation for different algorithms. Total energy consumption and number of orthogonal frequency channels 
needed to accomodate all the communication requests within the cell. Values are averaged over 1000 random configurations. 


VII. Conclusions 

We investigated the problem of energy efficient mode selection and resource allocation for 
network-assisted D2D communications in a dynamic TDD system. The problem has been anal¬ 
ysed under two frequency channel allocation strategies (with and without interference among 
D2D pairs) and with two objectives (total user energy and total system energy). For each 
configuration we derived the optimal solution to the corresponding MINLP formulation. For 
the full orthogonality scenario, we demonstrated how the optimal solution could be obtained 
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in polynomial (and sometimes even linear) time. When D2D pairs interfere with each other, 
on the other hand, we found it much harder to compute the optimal solution in an efficient 
manner. A customized branch-and-bound solver was therefore complemented by a more practical 
low-complexity heuristic. Through numerical simulations, we found that network-assisted D2D 
communications can yield significant energy savings, and that the heuristic algorithm could find 
near-optimal solutions while respecting practical implementation constraints. 
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